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ABSTRACT 
Descriptions  of  several  digital  computer  programs  for  electronic 
circuit  analyses  are  given.     A  listing  of  a  special-purpose  (Class  "C" 
amplifier  design)  program  is  given.     Comparisons  of  several  general 
purpose  programs  are  made.     An  extensive  description  of  the  CADLIB 
is  included.     The  discussion  of  the  CADLIB  points  out  numerous 
modifications  to  the  general  purpose  programs.     These  modifications 
include  such  things  as  the  addition  of  printer-plots,   more  accurate 
algorithms,    etc.     Also  many  of  the  problems  encountered  in  adapting 
the  various  programs  to  the  NPGS  machine  (IBM  360  model  67)  are 
enumerated.     Appendicies  are  included  to  provide  a  reference  for  coding 
the  members  of  the  CADLIB  for  which  no  instruction  manual  was  avail-  . 
able.     An  extensive  bibliography  on  applied  computer  aided  design  and 
the  job  control  language  details  of  how  the  CADLIB  was  created  are 
included.     In  addition  two  extensive  examples  using  the  CADLIB  are 
included. 
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CHAPTER  I 


INTRODUCTION 


In  circuit  design  a  most  time-consuming,    error-prone,   and 
frustrating  chore  is  the  repetitious  solving  of  network  equations  or 
the  repetitious  measurements  taken  from  a  breadboarded  circuit  as 
some  parameter(s)  is  varied.     However,    the  digital  computer  thrives 
on  simple-minded  repetitious  work  while  attaining  perfection  in  the 
writing  of  equations  and  near  perfection  in  their  solutions.     The 
accuracy  of  the  solutions  may  be  even  further  enhanced  by  the  use  of 
double  precision  calculations  which,   under  the  present  trend  of  com- 
puter design,   is  very  nearly  as  fast  as  single-precision  calculation, 
e.g.  ,    the  ADD  operation  in  the  IBM  360  requires  2.43  microseconds 
in  single-precision  and  2.45  microseconds  in  double-precision. 

Another  point  in  regards  to  accuracy  is  that  no  solution  is  any 
more  accurate  than  its  mathematical  model  irrespective  of  the  number 
of  significant  digits  carried  in  the  computation  process.     Modelling 
of  resistors,    coils,   and  capacitors  with  the  necessary  fringing  or 
stray  effects  presents  few  problems.     However,   an  accurate  model 
for  an  MOS-FET  or  integrated  circuit  is  prohibitively  complicated 
for  the  slide- rule  or  calculator  approach.     But  by  making  use  of  the 
speed  of  the  computer  virtually  any  model  can  be  effectively  used. 

In  the  following  discussions  a  basic  knowledge  of  the  Fortran 
programming  language  is  assumed.     There  are  also  numerous 


references  to  OS/  360  job  control  language  but  the  understanding  of 
these  references  is  not  necessary  to  make  use  of  this  report.     In 
particular  the  discussions  on  how  the  various  programs  were  placed 
on  a  disk  is  intended  primarily  for  guidance  to  one  who  is  interested 
in  modifying  or  adding  to  these  programs. 
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CHAPTER  II 
COMPUTER  AIDED  DESIGN  PROGRAMS 

INTRODUCTION 

There  are  several  hundred  digital  computer  programs  presently 
in  use  for  circuit  design.     These  programs  run  the  gamut  from  simple 
specialized  programs  for  a  particular  circuit  configuration  to  extremely- 
complicated  packages  that  have  widely  varying  capabilities  such  as: 
sensitivity  coefficients,   worst  case,   DC,   AC,   and  transient  analysis, 
plotted  outputs,   and  recently  on-line  graphics  (light  pen  and  CRT). 
Appendix  A  contains  an  example  of  a  specialized  program  to  "optimize" 
a  Class  "C"  amplifier  feeding  a  resistive  load.     As  is  characteristic 
of  these  specialized  programs  the  Class  "C"  program  is  of  limited 
value  and  the  input  is  very  restrictive  both  in  content  and  form.     Several 
worthwhile  specialized  programs  in  various  Fortran  languages  have 
been  listed  and  documented  in  JPL's  Technical  Report  32-9*79, 
references  18  and  19.     Examples  of  the  general  purpose  type  of 
program  comprize  the  remainder  of  this  report. 

GENERAL  DISCUSSION 

As  pointed  out  above  the  capabilities  of  the  available  programs 
vary  so  widely  in  applicability  and  scope  that  any  comparison  must  be 
limited  to  the  few  large  general-purpose  programs.  Table  1  lists  the 
major  characteristics  of  some  of  the  major  programs.     The  input  to 
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1  2 

these  programs  is  usually  engineer-oriented    and  free-format    ,   i.e.  , 

not  Fortran.     Typically  the  user  is  required  to: 

(1).     Number  the  nodes  sequentially  from  zero  or  one, 

(2).     Either  number  the  branches  sequentially  or  give  them  a 

unique  name,   and 

(3).     Code  the  branches  on  cards  using  commas  as  delimiters. 

The  limitations  of  these  programs  must  also  be  considered. 
All  state- variable  based  programs  suffer  from  two  common  limitations 
These  are: 

(1).     Numerical     errors  including  round-off  errors  and 

truncation  errors  ,   and 
(2).     The  time  constant  problem  generated  by  the  network 
being  analyzed. 

The  first  limitation  is  caused  by  the  limited  number  of  digits  in  a  com- 
puter word  and  can  be  minimized  by  the  use  of  better  algorithms  and 
double-precision  arithmetic.     The  second  limitation  is  the  result  of 
solving  continuous  integral  equations  by  quantized  numerical  methods. 
To  obtain  accuracy  of  integration  a  small  step  size  is  needed  for  even 

By  engineer-oriented  we  mean  that  a  branch  is  described  by  its 
node  numbers  and  a  positive  direction  of  current  flow. 

2 

By  free  format  we  mean  that  specific  columns  on  the  card  are 

meaningless  and  delimiters  such  as  commas  are  used  to  separate  the 
data. 
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ECAP 

NET-1 

CIRCUS 

SCEPTRE 

PREDICT 

NASAP 

Init.    Cond. 

X 

X 

X 

X 

X 

X 

Sens.    Coeff. 

X 

X 

X 

X 

DC  Analysis 

X 

X 

AC  Analysis 

X 

Trans.   Anal. 

X 

X 
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X 

X 

X 

Max  Element 

200 
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300 

300 
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50 
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300 

100 
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X 

X 
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xl,2 
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x1-2 
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x1'2 

xl,2,3 
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FTN 
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FTN 
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FTN/ FAP 

Deck  Size 
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30,000 
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15,000 

10,000 
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Nodal 

St  Var 

St  Var 

St  Var 

St  Var 

Dich 

Integration 

Imp. 

Pre. 

Exp. 

Exp,  Trap 

Trap 

Numer . 

Corr . 

&  R-K  4 

Notes: 

Ca 

n  be  consta 

.nt 

2 

Ca 

n  be  functi< 

on  of  time , 

table,    o 

r  equation 

3 

Ca 

n  be  a  func 

tion  of  a  si 

ibroutine 

Table  1.       A  comparison  of  the  Major  Programs. 
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the  smallest  time  constant  of  the  circuit.     But  to  obtain  the  complete 
transient  or  steady- state  solution  the  integration  must  be  carried  out 
to  about  five  of  the  largest  time  constants  of  the  circuit.     Thus  we  once 
again  hit  that  perennial  "trade- off"  problem  and  about  the  only  solution 
at  the  present  is  to  eliminate  as  many  of  the  small  time  constants  as 
practical  and/ or  use  a  variable  step- size  integration  technique  to 
minimize  the  error. 

Another  method  of  circumventing  the  first  limitation  is  to  use 
a  topological  approach.     However,   with  growing  circuit  size  computer 
run  time  rapidly  becomes  excessive  due  primarily  to  the  requirement 
that  all  possible  trees  must  be  found.     This  tree-finding  problem  will 
be  treated  more  thoroughly  in  a  later  section  dealing  with  the  CALAHAN 
program.     At  the  moment  topological  methods  are  limited  to  networks 
of  25-30  branches  and  15-20  nodes. 

ECAP  (Electronic  Circuit  Analysis  Program) 

ECAP    Cl6j    was  developed  through  a  joint  effort  of  IBM  and 
the  Norden  division  of  United  Aircraft  Corporation  and  released  in 
June  19^5  for  the  IBM  1620  computer.     It  has  since  been  translated  for 
use  on  the  IBM  7040/44,   7090/94,   and  the  System  360  series.     ECAP 
with  modifications  is  probably  the  most  widely  used  general  purpose 
CAD  (Computer  Aided  Design)  program  today.     It  is  one  of  very  few 
(if  not  the  only)  general  purpose  programs  to  offer  both  DC  and  AC 
analysis  with  an  automatic  parameter  modification  feature  in  both  its 
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DC  and  AC  analyses.     ECAP  will  provide,   as  output,   the  sensitivity 
coefficients,   worst  case,   and  standard  deviations  all  of  which  are 
noticeably  lacking  in  many  other  general  purpose  programs.     Another 
positive  feature  is  the  switch  (not  available  on  most  other  programs) 
which  provides  a  very  flexible  piece-wise  linear  modelling  capability 
to  approximate  non-linear  elements. 

A  major  limitation  of  ECAP  is  its  lack  of  a  non-linear  capability 
except  for  its  extremely  good  piece- wise  approximation  capability. 
The  major  problems  in  using  ECAP  are  its  lack  of  a  model  library  for 
frequently  used  components  and  the  problem  of  data  processing,    e.g.  , 
a  single  node  voltage  is  unobtainable  but  must  be  separated  from  a 
tabulated  list  of  all  node  voltages. 

NET-1 

NET-1      flOJ     was  developed  at  Los  Alamos  Scientific  Labora- 
tory of  the  University  of  California  at  Los  Alamos,   New  Mexico,   under 
the  auspices  of  the  United  States  Atomic  Energy  Commission.     NET-1, 
like  ECAP,   is  a  general  purpose  program  but  suffers  most  noticeably 
from  the  lack  of  an  AC  capability  and  the  unfortunate  (?)  choice  of 
FAP  as  a  primary  programming  language;  thereby,   limiting  its  users 
to  the  IBM  7090/94  and  Maniac  II  computers  as  conversions  to  other 
machines  would  be  a  major  undertaking.     However,   NET-1  does 
have  a  built-in  non-linear  model  for  transistors  and  diodes  but  these 
models  require  some  36  parameters  for  the  transistors  and  13  for 
the  diode  and  the  model  must  be  entered  before  it  can  be  used.     These 


15 


parameters  are  not  readily  available  from  manufacturers  and  must  be 
measured  from  several  devices  of  each  type  to  obtain  typical  values. 
NET-1  became  available  in  October  1962  and  has  since  been  rather 
widely  distributed  with  its  users  swapping  semiconductor  character- 
istics rather  freely. 

CIRCUS 

The  CIRCUS  program    ]J.Z  &   1  3 J  was  developed  at  the  Boeing 
Company,    Seattle,    Washington,   primarily  for  use  in  studying  the 
effects  of  nuclear  radiation  on  electronic  circuits,   but  its  applicability 
to  transient  analysis  of  all  types  led  to  its  use  as  a  general  purpose 
program.     CIRCUS  has  the  modeling  capability  (and  limitations)  of 
NET-1  but  was  written  in  Fortran  IV  (98%)  and  Assembly  (2%)  languages, 
and  the  conversion  to  the  System  360  has  been  made  by  Lincoln  Labora- 
tories at  MIT.     CIRCUS  is  relatively  fast  and  accurate  but  suffers 
from  the  same  fixed  models  as  NET-1.     This  can  lead  to  unnecessarily 
long  run  times  due  to  the  small  internal  time  constants  of  some  charge- 
control  models. 

SCEPTRE 


Like  CIRCUS  a  desire  for  a  transient  analysis  of  the  effects 
of  nuclear  radiation  led  to  the  development  of  SCEPTRE    [1 1 J    .     It 
was  developed  by  IBM  for  Kirkland  AFB.     One  of  the  contract  stipula- 
tions was  "100%  Fortran"  and  this  led  to  a  very  machine  dependent 
program  completely  defeating  the  purpose  of  the  Fortran  restriction. 
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SCEPTRE  has  two  unique  features  not  found  in  other  programs. 

The  first  is  the  ability  to  scan  the  loop  matrix  of  the  circuit, 
select  an  efficient  partition,    generate  a  Fortran  program  describing 
the  circuit,   and  then  execute  the  generated  program.     If  desired,   this 
Fortran  program  can  be  saved  and  executed  at  a  later  time  or  in  steps 
as  might  be  desired  when  the  run  time  becomes  long.     Parameter 
modifications  make  changes  only  in  this  generated  program  so  that  the 
load  time  (time  spent  generating  the  Fortran  program)  is  encountered 
only  once  for  each  circuit.     This  load  time  originally  required  4-5 
minutes  of  run  time;  however,   this  has  been  lowered  to  about  2  minutes 
[2\]  through  the  use  of  more  efficient  overlays. 

The  second  unique  feature  of  SCEPTRE  is  its  tremendous 
modelling  flexability.     SCEPTRE  will  accept  and  store  any  circuit  in 
its  model  library.     The  connections  and  their  order  is  specified  at 
the  time  of  entering  the  model  into  the  library  and  a  simple  call  will 
connect  the  model  anywhere  it  is  desired.     The  value  of  the  model 
library  for  "standard"  or  often  used  circuits  is  readily  apparent.     Due 
to  its   size  and  its  modelling  capability  SCEPTRE  is  probably  potentially 
the  most  powerful  general-purpose  CAD  program  available  today. 

PREDICT 

PREDICT  is  the  forerunner  of  SCEPTRE  but  as  originally  for- 
mulated it  has  no  stored  model  capability  and  had  a  rather  simple- 
minded  trapozoidal  integration  algorithm.     However,   an  IBM  360 
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version  has  been  written  and  some  of  the  features  of  SCEPTRE  were 
added  to  PREDICT  to  form  PREDICT- 360. 

NASAP  (Network  Analysis  for  System  Application  Program) 

NASAP  £4  &  7  j  is  mentioned  primarily  because  it  has  had  a 
maximization  of  advance  publicity  in  spite  of  the  fact  that  it  is  not 
operational  on  any  machine  at  the  moment.     But  NASAP  is  planned  to 
have  virtually  all  of  the  features  of  most  of  the  available  general-purpose 
programs  available  today.     NASAP  requires  no  lengthy  matrix  inversions, 
thereby  eliminating  the  primary  source  of  round- off  and  truncation 
errors.     NASAP  utilizes  a  dichotomous  approach  with  a  topological 

analysis  of  many  small  circuits  and  preliminary  testing  indicates  that 

5 
the  computer  run  time  is  reduced  by  factors  as  large  as  10     from  the 

state  variable  based  programs. 
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CHAPTER  III 


CADLIB 


The  CADLIB  is  a  computer  library  of  computer  programs 

designed  to  aid  the  designer  in  the  analysis  of  various  circuits.     The 

CADLIB  is  presently  maintained  in  executable  form  on  a  direct  access 

disk  named  J0BOO3  in  the  computer  center  of  the  Naval  Postgraduate 

School.     Support  for  the  CADLIB  is  furnished  in  the  form  of  a  sequential 

data  set  named  MODEL  (also  on  J0BOO3)  for  use  as  a  model  library  for 

one  of  the  members  of  CADLIB.     The  various  members  of  CADLIB  are 

executable  members  of  a  partitioned  data  set  and  they  consist  of  the 

following: 

ECAP  slightly  modified  by  J.  G.   Sketoe  to  furnish 

printer  plots      L^J   • 

CIRCUS  OS/  360  version  distributed  by  MIT  except 

that  the  default  model  unit  1 1  capability 
was  inserted. 

CALAHAN  slightly  modified  by  J.  G.  Sketoe  to  furnish 

automatic  printer  plots  both  of  freq  and 
transient  responses. 

POTTLE  slightly  modified  by  J  .   G.   Sketoe  to  furnish 

automatic  printer-plots  of  the  freq  response. 
Also  the  matrix  inversion  has  been  changed 
to  double-precision  arithmetic. 

POTTOL  the  POTTLE  program  with  a  free-field  format 

capability  as  designed  by  J.  G.   Sketoe 
(See  Appendix  D) . 

MODEL  the  CIRCUS  modeling  program  as  described  in 

reference  13.     Operates  on  data  set  MODEL. 
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The  CADLIB  data  set  was  constructed  using  a  SYSLM0D  of 

the  following  form: 

//SYSLM0D  DD  DSNAME=  CADLIB,  UNIT=  2  311,  DISP=(NEW, KEEP) ,     X* 

//  V<Z)LUME=SER=J0BOO3,  X 

//  SPACE=(CYL,(25,2,2)),  X 

//  LABEL=EXPDT=99350 

And  the  individual  members  were  added  using 

//SYSLM0D  DD  DSNAME=  CADLIB,  UNIT=2311,  DISP=0LD,  X 

V0LUME=SER=J0BOO3 

Execution  of  the  various  members  requires  different  sequences 

of  control  cards.     These  sequences  will  be  given  below  as  each  member 

is  discussed. 

ECAP 

The  first  program  obtained  for  the  CADLIB  was  IBM's  Elec- 
tronic Circuit  Analysis  Program.     The  360  version  requires  approx- 
imately 290K  bytes  of  core  when  used  without  overlays  as  in  CADLIB. 
ECAP  was  obtained  from  Share,   a  facility  of  IBM's,   in  the  E-level 
compiler  edition.     As  the  OS/ 360  version  was  ordered  only  the  360 
manual  [20J  was  included.     However,   this  manual  only  pointed  out  the 
differences  in  ECAP- 1620  and  ECAP- 360.     A  manual  was  then  ordered 
for  the  1620  version  Q.oj  and  two  months  later  it  arrived.     In  the  mean- 
time by  using  ECAP  problems  published  in  various  journals  it  was 


This  X  as  in  all  of  the  following  control  cards  must  be  in  column  72. 
Also  the  first  character  of  the  following  card  (except  for  the  double 
slashes)  must  be  in  column  16. 
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discovered  that  the  received  version  simply  would  not  work  properly. 

The  major  problems  and  their  solutions  were: 

(1).     The  release  9  G-level  "in-house"  compiler  failed  to  con- 
struct the  proper  instructions  for  overlay.     At  that  time 
the  computer  center  policy  was  to  run  two-256K  byte 
machines  rather  than  one  512K  byte  machine  so  that 
various  hardware  testing  could  be  done  concurrently  with 
computer  service.     Use  of  the  older  but  slower  H-level 
compiler  cured  this  trouble.     (This  fault  of  the  G- compiler 
has  been  cured  in  later  releases). 
(2).     The  next  problem  discovered  was  that  the  language 

processor  portion  of  ECAP  did  not  function  correctly 
when  compiled  on  either  the  G-  or  H-level  compiler.     The 
cause  of  this  trouble  lay  at  the  fundamental  differences 
in  the  E-level  Fortran  compiler  and  the  G/H  level  com- 
pilers.    As  the  language  processor  converts  an  engineer- 
oriented  circuit  description  into  the  required  computer 
description,   literal  data  (alphanumeric  characters)  is 
manipulated  vice  purely  binary  numbers  (floating  point 
numbers  are  a  special  coding  of  binary).     This  manipula- 
tion is  done  by  using  a  unique  combination  of  bits  to 
represent  a  particular  character.     When  an  arithmetic  IF 
statement  is  used  to  decode  this  data,   the  data  is  treated 
as  numbers  and  not  as  literal  characters.     For  example, 
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the  Fortran  statement 

IF  (M-N)        10,    15,    10 
compiles  into  an  arithmetic  subtraction  of  N  from  M, 
examine  the  resultant,   then  branch  to  statement  10,    15, 
or  10  according  to  whether  the  resultant  is  negative,    zero, 
or  positive.     The  contents  of  a  computer  word  or  cell  is 
considered  positive  or  negative  according  to  whether  the 
first  bit  is  a  zero  or  one,    respectfully.     In  the  IBM  7090/94 
machines  the  literal  data  codes  were  mixed  as  positive 
and  negative  numbers  with  the  special  characters  (,  .  etc.) 
intermingled  within.     This  mixture  was  observed  in  the 
E-level  machines  while  in  later  machines  this  coding  was 
changed  such  that  alphanumerics  are  negative  (and 
sequential)  and  special  characters  are  positive.     The 
solution  to  this  problem  was  quite  obvious  once  the  source 
of  the  error  was  known.     All  literal  data  arithmetic  IF's 
were  replaced  with  logical  IF's  which  do  their  comparisons 
on  a  bit-for-bit  basis  rather  than  an  arithmetic  answer. 
(3).     Solving  of  these  two  problems  finally  got  our  copy  of  ECAP 
at  least  partially  working.     The  next  problem  encountered 
was  that  the  automatic  parameter  variation  feature  did 
not  function  properly.     A  cry  for  help  to  IBM  was  answered 
with  a  list  of  errors  and  omissions  in  the  originally  dis- 
tributed version  of  ECAP.     The  appropriate  entries  cured 

this  problem . 
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(4).     Next  an  addressing  error  was  encountered  when  trying 
to  obtain  a  frequency  response  solution  with  an  overlaid 
program.     As  the  various  control  sections  in  overlays 
complicate  debugging,    an  attempt  was  made  to  find  the 
error  while  using  all  512K  bytes  of  core  without  overlays. 
The  addressing  error  disappeared.     Reinsertion  of  the 
overlay  brought  the  bug  back.     The  only  conclusion  was 
that  the  bug  was  in  the  machine  constructed  overlay. 
Next  run  time  comparisons  were  made  on  overlay  vs. 
non-overlay  versions.     The  results  were  that  overlaid 
programs  required  2-4  times  as  much  run  time  as  non- 
overlay.     By  this  time  all  51 2K  was  being  used  for  almost 
all  runs  and  the  decision  to  use  ECAP  without  overlay  in 
the  CADLIB  was  made.     The  bug  was  brought  to  the  attention 
of  the  computer  center  personnel  and  left  there. 

Presently  all  known  bugs  have  been  worked  out  of  the  NPGS 
version  of  ECAP  in  addition  to  the  addition  of  a  printer  plotting  addition 
as  explained  in  reference  17.     As  User's  Manuals  are  in  short  supply 
an  abbreviated  manual  has  been  written  and  is  published  as  reference  17. 

The  required  control  cards  to  execute  the  CADLIB  version  of 
ECAP  are  shown  in  figure  1  . 
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//ECAP      JOB    (EAC6.01FT),1     NAME     '  ,MSGLEVEL=1 
//J0BLIB  DD  DSNAME=CADLIB,UNIT=2311  ,DISP=OLD,  X 

//  VOLUME=SER=J0BOO3 

//  EXEC      PGM    ECAP 

//FTO6F001  DD      SYS0UT=A 

//FTO5F001  DD      * 

Batched  ECAP  coded  input  decks 
/* 


Figure  1.       Control  Cards  to  Execute  ECAP. 

Basically  ECAP  is  formulated  to  fill  the  nodal  conductance 
matrix  as  the  input  branches  (B-cards)  are  read  in.     The  various 
command  cards  then  set  the  appropriate  flags  to  indicate  which 
algorithms  are  executed  to  obtain  the  desired  outputs.     This  conductance 
matrix  is  available  as  an  output  request  under  the  MISCELLANEOUS 
request. 

One  additional  note  worthy  of  comment  about  ECAP  is  that  it 
only  interprets  the  first  two  characters  of  any  word,    e.g.  ,   PR  is  the 
same  as  PRINT,   V0    is  the  same  as  V0LTAGES,   BE  is  the  same  as 
BETA,    etc. 

CIRCUS 

A  copy  of  the  CIRCUS- 360  program  was  obtained  from  Lincoln 
Laboratory,   MIT,    Lexington,   Massachusetts.     As  with  ECAP  the 
received  version  would  not  work  although  MIT  had  exactly  the  same 
type  machine  as  NPGS,   i.e.  ,   a  model  65/67  IBM  360.     The  MIT  des- 
cription of  the  distributed  CIRCUS  tape  is  given  in  Appendix  B.     The 
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various  problems  encountered  with  our  copy  were: 

(1).     The  CIRCUS  tape,   as  received,   had  an  H-level  compiler 
output  version  on  file  #  1  (Appendix  B)  which  should  have 
loaded  directly  from  the  tape.     However,   a  tape  error 
had  occurred     and  the  tape  had  to  be  recompiled. 

(2).     At  MIT  there  exists  three  resident  routines  that  access 
the  day,   year,   hour,   minute,   and  second.     These 
routines  were  used  by  CIRCUS- 360.     As  these  routines 
did  not  exist  at  NPGS  the  initial  solution  was  to  provide 
dummy  routines  to  get  CIRCUS  operating.     While  this 
solution  worked  it  eliminated  the  variable  time  limit  of 
CIRCUS  as  explained  in  the  User's  Manual     Q2J     .     The 
author  then  proceeded  to  learn  enough  assembly  language 
to  write  a  set  of  routines  to  accomplish  this.     These 
routines  are  now  added  to  our  version  of  CIRCUS  as  well 
as  entered  in  the  resident  library  of  our  computer  center, 

(3).     While  translating  CIRCUS  MIT  evidently  inserted  a  de- 
bugging print- out  at  each  integration  step.     As  this  led 
to  excessive  and  useless  print-out  of  the  integration  step 
time  and  size,   both  in  time  and  paper  this  routine  was 
slightly  modified  to  eliminate  this  waste. 


A  bEND  card  had  become  BEND  where  "b"  implies  a  12-2-9  punch, 
A  "B"  character  implies  a   12-2  punch. 
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(4).     No  overlay  table  was  included  in  CIRCUS- 360  as  it  will 
fit  into  a  256K  byte  machine  if  the  load  point  is  below 
B000  (See  Appendix  B) .     However,   the  load  point  at  NPGS 
is  B820    ,  with  just  the  "in-house"  G-compiler  thereby 
necessitating  an  overlay.     This  overlay  was  constructed 
and  is  shown  in  figure  2.     This  overlay  was  used  but  was 
soon  dropped  for  the  following  reasons: 

(a).    Substantial  decrease  in  run  time. 

(b).    Prior  experience  with  ECAP  indicated 

errors  might  occur, 
(c).    Using  the  256K  byte  machine  only  158K 
was  available  with  HASP.     CIRCUS  even 
with  this  overlay  required  153K  leaving 
almost  no  room  for  increased  compiler 
size  (which  came  later), 
(d).    By  going  to  512K  we  were  able  to  increase 
our  maximum  circuit  size  to  circuits  of 
over  70  equations  (Appendix  B). 

(5).     The  MIT  description  lists  two  known  bugs  (Appendix  B) 

in  their  program.     The  second  one  of  these  failed  to  show 
up  in  our  recompiled  version  of  CIRCUS.     This  led  to  a 
successful  attempt  to  correct  the  tape  error  mentioned 
earlier  so  that  this  bug  could  be  checked  out.     As  indicated 
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Figure  2.       CIRCUS- 360  Overlay  Requiring 
approximately  153K  bytes. 
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in  Appendix  B  the  MIT  version  failed  to  recognize  a 
card  of  the  form 

USE  MODEL  LIBRARY,   n 
where  n  is  integer  denoting  the  logical  unit  of  the  model 
library.     The  NPGS  version  did  recognize  the  card 
correctly.     Naturally  the  NPGS  version  was  placed  in 
CADLIB . 
(6).     An  error  of  omission  by  MIT   but   not  mentioned  by  them 
in  their  write-up  was  that  the  default  to  unit  1 1  (page  25  of 
User's  Guide  Q2j  )  was  omitted.     This  has  been  corrected 
in  the  NPGS  version. 

Presently  CIRCUS- 360  is  in  CADLIB  and  is  working  properly. 
As  CIRCUS  uses  several  different  logical  units  a  cataloged  procedure 
has  been  established  for  CIRCUS  in  order    to  simplify  the  job  control 
card  requirements.     The  standard  execution  of  CIRCUS  using  logical 
unit  11  for  the  model  library,   when  needed,   is  shown  in  figure  3. 


//CIRCUS      J0B      (EAC6,01FT),«    NAME    ',MSGLEVEL=1 
//J0BLIB  DD  DSNAME=  C ADLIB,  UNIT=2 311  ,DISP=<Z>LD,  X 

//  V0LUME=SER=J0BOO3 

//        EXEC         CIRCUS 
//GO.SYSIN      DD      * 

Batched  CIRCUS  coded  input  decks 
/* 


Figure  3.       Standard  CIRCUS  Job  Control  Language 
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If  logical  assignments  other  than  or  in  addition  to  the  default 
assignments  are  desired  (such  as  use  of  "SAVE  TAPE"),   these  units 
must  be  defined  in  a  manner  similar  to  that  shown  in  figure  4  where 
unit  15  is  defined  for  use  as  a  "SAVE  TAPE'  or  as  a  model  library  tape 

A  listing  of  the  actual  cataloged  procedure  to  execute  CIRCUS 
is  given  in  figure  5  for  reference  althouth  it  should  not  be  needed 
except  to  perhaps  extend  or  replace  it. 


//CIRCUS     J0B  (EAC6,01FT),'  NAME 

1 ,  MSGLEVEL=  1 

//J0BLIB  DD  DSNAME=CADLIB,UNIT= 

2311,DISP=0LD. 

X 

// 

V0LUME=SER=J0BOO3 

// 

EXEC 

CIRCUS 

IIG0 

.FT15F001     DD    DSNAME=  SAVE,  UNIT=  2400  ,  DISP= 

(NEW, KEEP), 

X 

II 

V0LUME=SER=XXX, 

LABEL=YYY, 

X 

II 

DCB=(BLKSIZE=  1 204 

,RECFM=V) 

IIGQ 

.SYSIN 

DD       * 

Batched  CIRCUS  coded  input  deck 

s 

/* 

Figure  4.     Defining  Logical  Unit  15  for  use  by  CIRCUS 


MEMBER  NAME       CIRCUS 

//GO                     EXEC       PGM=CIRCUS 

//FTO1F001    DD  DSNAME=SYS1.SYSUT1,DISP=0LD, 

X 

//                                      DCB=(BLKSIZE=1204,RECFM=V 

//FTO3F001    DD  DSNAME=SYS1.  SYSUT2  ,  DISP=0LD, 

X 

//                                     DCB  =  (BLKSIZE  1204,  RECFM=  V) 

//FT12F001    DD   DSNAME=SYS1.  SYSUT3,  DISP=0LD, 

X 

//                                      DCB=(BLKSIZE=1204,RECFM=V) 

//FT11F001    DD    DSNAME=M0DEL,UNIT=2311,DISP=0LD, 

X 

//                                      V0LUME=SER=J0BOO3 

//FT06F001  DD    SYS0UT=A 

//FT05F001  DD    DDNAME=SYSIN 

Figure  5.     Cataloged  Procedure  "CIRCUS" 
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The  formulation  of  CIRCUS  is  based  in  the  state  varible  approach 
as  advanced  by  Bryant     f 22J     and  extended  by  Wilson  and  Massena     ("23*] 
A  detailed,   although  not  rigorous,    explanation  of  the  mathematical 
theory  of  CIRCUS  is  given  in  the  Program  Manual    [^13*J     . 

MODEL 

The  MODEL  member  (alias  UPDATE)  of  CADLIB  is  a  program 
to  list,   modify,    generate,    etc.   the  model  library  as  described  in 
Chapter  4  of  the  CIRCUS  Program  Manual     [l3^J     .     The  CIRCUS  tape 
as  received  from  MIT  had  a  compiled  version  of  the  MODEL  program 
as  file  #2  (see  Appendix  B) .     This  file  was  used  as  input  to  the  linkage 
editor  to  form  member  MODEL  which  operates  on  the  actual  CIRCUS 
model  library.     The  control  cards  required  for  execution  of  the  member 
MODEL  is  shown  in  figure  6.     The  model  library  is  described  by  the 
FT11F001   card  in  figure  6  which,   as  shown,   describes  a  date-protected 
sequential  data  set  named  MODEL  which  resides  on  the  same  JOB003 
disk  as  the  CADLIB.     If  one  wishes  to  use  the  program  or  member 
MODEL  to  operate  on  some  other  data  set  (i.e.  ,   magnetic  tape)  only 
this  FT11F001  must  be  changed. 
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//MODEL      JOB    (EAC6.01FT),  '  NAME  ',MSGLEVEL=1 

//J0BLIB  DD  DSNAME=CADLIB,UNIT=2  311  ,DISP=0LD, 

//  V<2>LUME=SER=J0BOO3 

//  EXEC        PGM=M0DEL  or         UPDATE 

//FT11F001      DD      DSNAME=M0DEL,UNIT=2311,DISP=0LD, 

//  V0LUME=SER=JOBOO3 

//FT12F001      DD      DSNAME=SYS1  .  SYSUT1 ,  DISP=0LD 

//FT06F001      DD      SYSOUT=A 

//FT05F001      DD      * 


X 


X 


Model  input  as  per  Chapter  4,    CIRCUS 
Program  Manual    £l  3J 


/* 


Figure  6.       Use  of  the  Model  Program. 


POTTLE 

The  copy  of  C.   Pottle's  Linear  Network  Analysis  Program  as 
received  required  only  changes  in  the  logical  unit  assignments  to  make 
it  compatible  with  the  NPGS  machine.     To  make  the  use  of  this  program 
feasible  from  the  CADLIB  all  outputs  were  assigned  to  unit  #  6  (printer) 
and  the  input  was  assigned  to  unit  #5  (card  reader).     However,   flexibility 
was  added  to  the  program  by  making  the  unit  designations  variables 
in  the  individual  read  and  write  statements  and  then  assigning  these 
variables  in  a  central  location  in  the  main  program.     These  variables 
are  then  transferred  between  various  subroutines  via  labeled  common. 
This  provides  the  user  with  the  facility  of  writing  out  some  data 
separately  from  other  data.     The  variable  names  and  their  assignments 
are  shown  in  figure  7. 
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Name 

Unit 

Use 

INTAPE 

5 

User's  input  data 

IOUTP 

6 

Copy  of  data  set  &  errors 

STATP 

6 

State  equations 

TRANTP 

6 

Transfer  functions,   poles, 

&   zeros 

FREQTP 

6 

Frequency  response 

Figure  7.       Default  POTTLE  I/O  Assignments. 

The  control  cards  required  to  execute  the  POTTLE  program 
from  the  CADLIB  is  shown  in  figure  8. 


//P0TTLE      J0B    (EAC6,01FT),M  NAME  *,MSGLEVEL=1 
//J0BLIB  DD    DSNAME=CADLIB,UNIT=2311,DISP=0LD, 
//  V0LUME=SER=JOBOO3 

//        EXEC  PGM=P0TTLE 

//FTO6F001        DD      SYS0UT=A 
//FT05F001         DD      * 


X 


/* 


POTTLE  coded  input    (See  Appendix  C) 


Figure  8.       Control  Cards  to  Execute  POTTLE. 

The  procedure     £l43   implemented  by  Mr.    Pottle  is  based  on  a 
recursive  algorithm  developed  by  Derisoglu    |^6j     which  extends  Bryant's 
algorithm  to  nonlinear  active  elements.     The  state  equations  are  formula- 
ted; then  the  frequency  response  is  determined  via  the  Souriau-Frame 
algorithm  to  obtain  the  matrix  (sI-A)    .       The  time  responses  are 
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obtained  by  factoring  the  characteristic  polynomial  and  performing 
partial  fraction  expansions  of  each  element  of  the  transfer  function 
matrix  or  by  direct  integration  of  the  set  of  first-order  equations. 

The  instructions  for  coding  the  POTTLE  program  are  included 
as  Appendix  C. 

Changes  to  the  POTTLE  program  were  limited  to  (1)  providing 
extended  accuracy  through  the  use  of  a  double-precision  vice  the 
original  single-precision  matrix  inversion  routine  and  (2),   the  addition 
of  an  automatic  plot  of  the  time  response  and  the  linear  frequency 
response. 

POTTOL 


The  frustrations  encountered  in  using  the  POTTLE  program 
pointed  out  the  desirability  of  a  translator  so  that  a  free-format  input 
could  be  used.     The  author  has  designed  such  a  translator  and  combined 
it  with  the  POTTLE  program  described  above  to  form  the  POTTOL 
program.     While  developing  the  Translator  several  other  convenience 
features  were  included. 

While  the  use  of  the  POTTOL  input  format  is  relatively  simple, 
a  complete  description  of  the  Translator  is  not  so  simple.     The  complete 
description  of  the  Translator  input  coding  is  given  in  Appendix  D.     The 
control  cards  to  implement  POTTOL  are  shown  in  figure  9  and  are 
very  similar  to  those  used  for  ECAP  and  POTTLE. 
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//P0TT0L      J0B    (EAC6.01FT),1  NAME  ',MSGLEVEL=1 
//J0BLIB    DD    DSNAME=CADLIB,UNIT=2311,DISP=0LD,  X 

//  V0LUME=SER=J0BOO3 

//  EXEC  PGM=P0TT(Z>L 

//FT06F001      DD  SYS0UT=A 

//FT05F001      DD  * 

Batched  POTTOL  coded  input  (See  Appendix  D) 
/* 


Figure  9.       Control  Cards  to  Execute  POTTOL. 


CALAHAN 


Probably  the  most  accurate  program  in  the  CADLIB  is  the 
CALAHAN    [33     which  was  developed  by  D.  A.   Calahan  while  at  the 
University  of  Illinois.     As  received,   the  program  was  coded  in  IBM 
7090  Fortran  II.     The  problems  encountered  with  the  program  were 
generated  primarily  by  a  complete  lack  of  knowledge  of  the  intricacies 
of  Fortran  II.     But  the  differences  were  not  as  great  as  feared.     The 
names  of  the  Fortran  library  routines  were  slightly  different  (e.g.  , 
COSF  vs.    COS)  and  the  passing  of  arguments  through  subroutines  could 
be  handled  differently  causing  the  computer  to  get  lost  when  passed 
in  the  same  manner  in  Fortran  IV.     The  major  problem  encountered 
in  implementing  and  testing  the  program  was  the  intepretation  of  the 
instruction  manual.     Appendix  E  is  a  set  of  instructions  which  (hopefully) 
will  clearify  the  coding  of  the  input  to  the  CALAHAN  program. 

The  CALAHAN  program  is  a  topologically  based  program  which 
is  the  primary  source  of  its  accuracy.     This  approach  evaluates  the 


34 


determinants  and  cofactors  of  the  node  admittance  matrix  by  finding 
the  related  n-tree  products  of  the  corresponding  linear  graph.     Basically 
the  matrix  solution  for  the  node  voltages  is  transformed  by  the  Binet- 
Cauchy  theorem    |  15  J   and  the  Maxwell's  formula  to  a  problem  of 
determining  all  tree  and  associated  n-tree  products.      This  tree-finding 
problem  quickly  replaces  the  integration  time  constant  problem  of  the 
state- variable  program  as  a  computer  run  time  limitation.     Regardless 
of  the  tree-finding  scheme  used  the  possible  combinations  of  branches 
increase  factorially  (about  the  most  rapidly  increasing  function  known) 
and  the  number  of  trees  increase  exponentially  (about  the  second  most 
rapidly  increasing  function  known).     For  example,    the  circuit  of  figure 
10  using  the  hybrid-pi  transistor  model  results  in  a  circuit  containing 
eleven  nodes  and  twenty- eight  branches  which  forms  456,  172  trees! 
While  all  topologically  based  programs  have  this  tree-finding  limitation, 
they  do  not  suffer  from  the  round- off  error  so  prevalent  in  the  state- 
variable  formulations. 
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Figure  10.       Wide-band  R-C  Coupled  Amplifier 
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The  execution  of  the  CALAHAN  program  from  the  CADLIB  is 
shown  in  figure  11  and  again  the  similarity  to  the  execution  of  the 
ECAP  and  POTTLE  should  be  noted. 


//CALAHAN      JOB    (EAC6,  01FT,  '  NAME  ' ,  MSGLEVEL=  1 
//J0BLIB  DD    DSNAME=  CADLIB, UNI T=2 3 11,DISP=0LD,  X 

//  V0LUME~SER=J0BOO3 

//  EXEC        PGM=  CALAHAN 

//FT06F001      DD      SYSOUT=A 
//FT05F001      DD      * 


/* 


Batched  CALAHAN  coded  input    (See  Appendix  E) 


are: 


Figure  11.       Control  Cards  to  Execute  CALAHAN. 
The  extensions  and  additions  made  to  the  CALAHAN  program 

(1)  An  automatic  printer-plot  of  the  requested  output  for 
both  frequency  response- calculations  and  time  response 
calculations  . 

(2)  A  double-precision  root-finding  algorithm  which 
usually  works  better  than  the  original  "Method  of 
Mueller"  algorithm  has  been  added.     The  roots  are 
found  by  both  methods,   two  polynomials  are  then  con- 
structed from  the  two  sets  of  roots,   the  magnitude  of 
the  coefficients  are  compared  to  the  original.     The  set 
of  coefficients  with  the  least  error  determines  the 
selected  set  of  roots. 
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Where 

V. 
Z(s)  =     F(s)     =       y^- 
in 

And 

V 
T(s)  =    F(s)    =       — 

in 


Considering  now  the  case  of  non-ideal  inductors  and  capacitors 
we  can  define  a  new  function  or  "s"  in  which  the  impedances/ admittances 
become  (r      +  Ls)  /   (g      +  Cs)  where  r      and  g      are  the  dissipative 

L.4  Lv  J_;  \*> 

elements  of  the  inductors  and  capacitors,    respectively.     Now  defining 


^l   "      L  \    '   ~c 


We  have 


24ft)=.tCs.t.O         and       YcCs)  =  c(stScN) 

At  this  point  we  could  depart  from  the  more  rigoreous  approach 
and  make  the  assumption  that      $    ~   *&    ~    0     •  thereby,   making  life 
very  simple  by  drawing  the  conclusion  that  our  function  with  non-ideal 
elements  is  formed  by  the  substitution  of  (s  +  £    )  for  s  in  the  original 
function.     In  order  to  test  the  validity  of  such  an  assumption  the 
dissipation  of  capacitors  and  coils  were  measured  using  a  General 
Radio  Type   1650  impedance  bridge.     These  components  were  chosen 
at  random  with  the  only  criteria  being  that  a  "usage  range"  be  obtained. 
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The  results  of  these  measurements  are  tabulated  along  with  calcula- 
tions of     St     and     $     in  table  2.     As  is  apparent  from  the  table  the 


Type 

Nomial 

Measured      D( 

)     Cal. 

Comments 

Cap 

Electro 
lytic 

4  ufd 

4  ufd 

270 

1.7 

Ignored 
see  text 

Oil 

2  ufd 

1.93 

3.9 

0.0245 

Paper 

.0001 

.00093 

1.8 

0.0113 

Minimum 

Mica 

.0047 

.00446 

4.7 

0.0295 

Mica 

200  pf 

182  pf 

4.0 

0.0251 

Cermic 

250  pf 

258  pf 

1.0 

0.0628 

Maximum 

Coils 

2-        r-f 

700  uh 

710  uh 

280 

1.76 

1-        r-f 

200  uh 

206  uh 

265 

1.66 

Minimum 

RF- Choke 

58     uh 

105 

6.60 

Air -  Dux* 



19.5uh 

.36 

2.26 

Air -Dux"' 

742  uh 

1.2 

7.55 

Maximum 

Locally 

25.2uh 

.52 

3.25 

wound 

*  Air  -Dux 

is  a  trade 

name  of 

a  type  of 

coil  used  f 

or 

transmitt 

er-to-ante 

nna  lead 

coupling 

Table  2.       Coil  and  Capacitor  Measurements 

numerical  value  of     £      or     ^      did  not  vary  excessively  over  the     "usage 
range"  except  for  the  single  4  ufd  electrolytic  capacitor.     As  electrolytic 
capacitors  have  relatively  little  use  where  the  pole  and  zero  positions 
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are  critical,   this  particular  capacitor  was  subsequently  ignored. 
However,   the  approximation  that    %  ~£     must  certainly  be  held  invalid 
or  at  least  suspect  as  these  parameters  differ  by  factors  of  100. 

By  utilizing  the  digital  computer,   the  analysis  of  a  particular 
circuit  can  be  reduced  to  the  coding  of  a  few  cards.     Utilization  of 
the  CALAHAN  digital  computer  program  in  accordance  with  the 
instructions  in  Appendix  E  provides  the  desired  quantities  with  a 
minimum  of  druggery. 

The  two  circuits  used  for  this  problem  were: 

(1)  The  Pi-Section  of  figure  12  which  is  often  used  to  couple 
a  transmitter  to  an  antenna,   and 

(2)  The  Chebyshev  filter  of  figure  13. 
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Figure  12.       Pi- Section  Coupler 
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Figure  13.         Chebyshev  Filter 
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rc 

ohms 

V 
out 

V. 
in 

T(s)     =      V. 
in 

Z(s)    =       I. 
in 

.01 
.05 
.  10 
.50 

100 
20 
10 
2.0 

-1.01 ,   -.505  +  jl.321 
-1.05,   -.525  +  jl.313 
-1. 10,    -.550  +  jl.303 
-1.50,    -.750  +  jl .200 

-.5812,   -.2194  +  j 1.3065 
-.6265,   -.2368  +  jl. 3038 
-.6830,   -.2585  +  jl. 3000 
-1.132,   -.4342  +  jl .2565 

Table  3.     Pi-Section  Pole  Movement  (i-        =     0,   vary    £       ) 

L  7     °C 


Since  it  was  suspected  that  the  capacitors  would  be  substantially 
"more  ideal"  than  the  inductors  in  their  effects  on  the  circuit  the  first 
computations  were  to  vary     Q-       over  the  "typical  range"  and  determine 
the  functions'  pole  positions  of  the  Pi-Section.     These  are  shown  in 
Table  3.     Next  similar  computations  for  the  "typical"     £      range 
were  made.     These  results  are  shown  in  Table  4. 


k* 

rL 

T(s) 

V 
out 

V. 
in 

L              rL 

1.0 

-1.0, 

-1.0    + 

j   1.4142 

5.0 

-1.0, 

-1.5858, 

-4.414 

10.0 

-1.0, 

-1.228, 

-9.772 

Table  4.         Pi-Section  Poles  vs. 

As  can  be  seen  from  these  calculations  the  effects  of  the  typical 
inductor  on  the  movement  of  the  poles  is  some  100  times  greater  than 
the  typical  capacitor.     This  correlates  nicely  with  the  measured      £   S 


41 


The  computations  for  varying     £>       in  the  Chebyshev  filter 
of  figure  13  were  carried  out  but  the  effects  were  again  on  the  order 
of  2%  (See  table  5). 


&o 

Zeros 

Poles 

0.01 
0.05 
0.10 

-0.005    +      jl.514 
-0.025    +      jl .5136 
-0.05      +      jl.513 

-0.9819 

-0.2237      +      jl.140 

-1.217        +      jl.425 

-1.041 

-0.2446      +      jl . 1378 

-1.2266      +      jl.4032 

-1.116 

-0.2711      +      jl.1356 

-1.2375      +      jl.3757 

A 


*n 


O 


Table  5.       Chebyshev  Filter  (Poles  and  Zeros  vs 
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Figure  14.       Pi- Section  with  Non-Ideal  Elements. 

The  transfer  function  of  figure   14  was  calculated  and  the  poles 
determined  for  various  values  of  RC1  and  RC2.     The  particular  values 
were  10 ,    100  ,   and  1000  ohms  which  corresponded  to  a     of      of  0 .  1  , 
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0.01,   and  0.001,    respectfully;  thereby,   bracketing  the  measured 

typical  values.     With  each  of  these  values  RL  was  varied  over  a  very 

-7 
wide  range  of  10        to  100  and  pole  positions  calculated.     These  results 

are  tabulated  in  table  6  and  graphed  in  figure  15. 

As  indicated  by  the  graphs  the  pole  positions  certainly  do  not 
move  straight  in  the  negative  real  direction.     In  fact,   there  is  quite 
a  severe  vertical  movement.     However,   notice  that  the  poles  do  come 
back  down  to  level  with  the  ideal  at  RL  ~  2  which  is  in  the  midrange 
of  the  measured    &      ;  thereby,   in  some  sense  justifying  the  "direct 
left  shift"  as  a  ball  park  design  criteria. 

The  effects  of  the  capacitor  admittance  is  to  move  the  pole 
position  slightly  to  the  left  and   downward  while  the  inductive  resistance 
moved  the  pole  up  until  RL        L  (numerically)  and  then  the  effect  be- 
comes quite  the  opposite.     The  major  point  to  be  gleaned  from  this 
problem  is  that  the  usually  neglected  inductor  resistance  (in  filter 
design)  can  be  of  such  a  magnitude  to  cloud  the  analysis  (with  ideal 
elements)  in  doubt.     The  CADLIB  provides  a  tool  to  perform  the 
mathematical  druggery  so  that  these  parasitic  elements  can  be  included 
in  the  original  design. 

ON  COMPENSATIONS  IN  INTEGRATED  CIRCUIT  DESIGN 

Integrated  circuits  has  become  one  of  the  most  important  areas 
in  space  and  military  electronics.     As  a  consequence,   ingegrated 
circuit  design  has  become  one  of  the  major  branches  of  network 
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RL 

RC1= 

=RC2=10 

RC1=RC 2=100 

RC1=RC 2=1000 

lO"7 

-.5500 

+  j.8352 

-.5050 

+  j.8631 

-.5005 

+  j.8657 

10"6 

-.5500 

+  j.8352 

-.5050 

+  j.8631 

-.5005 

+  j.8657 

lO"5 

-.5500 

+  j.8352 

-.5050 

+  j.8631 

-.5005 

+  j.8657 

0.1 

-.6000 

+  j.8660 

-.5550 

+  j.8905 

-.5505 

+  j.8928 

0.2 

-.6500 

+  j.8930 

-.6050 

+  j.9143 

-.6005 

+  j.9163 

0.4 

-.7500 

+  j.9368 

-.7050 

+  j.9524 

-.7005 

+  j  .9558 

0.5 

-.8000 

+  j.9539 

-.7550 

+  j.9669 

-.7505 

+  j.9681 

0.7 

-.9000 

+  j.9798 

-.8550 

+  j.9879 

-.8505 

+  j.9886 

0.9 

-1.000 

+  j.9950 

-.9550 

+  j.9985 

-.9505 

+  j.9987 

1.0 

-1.050 

+  j.9987 

-1.005 

+  jl.000 

-1.000 

+  jl.000 

1.1 

-1.100 

+  jl.000 

-1.055 

+  j.9990 

-1.050 

+  j.9988 

1.2 

-1.150 

+  j.9987 

-1.105 

+  j  .9955 

-1.100 

+  j.9950 

1.3 

-1.200 

+  j.9950 

-1.155 

+  j.9895 

-1.150 

+  j.9888 

1.4 

-1.250 

+  j.9887 

-1.205 

+  j.9808 

-1.200 

+  j.9799 

1.5 

-1.300 

+  j.9798 

-1.255 

+  j.9695 

-1.250 

+  j.9684 

1.6 

-1.350 

+  j.9682 

-1.305 

+  j .9555 

-1.300 

+  j.9541 

1.7 

-1.400 

+  j.9539 

-1.355 

+  j.9386* 

-1.350 

+  j.9369 

1.8 

-1.450 

+  j.9368 

-1.405 

+  j.9187 

-1.400 

+  j.9167 

1.9 

-1.500 

+  j.9165 

-1.455 

+  j .8955 

-1.450 

+  j .8933 

2.0 

-1.550 

+  j.8930 

-1.505 

+  j.8689 

-1.500 

+  j.8663 

2.2 

-1.650 

+  j.8352 

-1.605 

+  j.8037 

-1.600 

+  j.8004 

2.4 

-1.750 

+  j.7599 

-1.705 

+  j.7190 

-1.700 

+  j.7146 

2.6 

-1.850 

+  j.6614 

-1.805 

+  j.6066 

-1.800 

+  j.6007 

2.8 

-1.950 

+  j.5268 

-1.905 

+  j.4461 

-1.900 

+  j.4369 

3.0 

-2.050 

+  j.3123 

-2.005 

+  j.0999 

-2.000 

+  j.0316 

3.2 

-2.470 

,   -1.830 

-2,555 

;  -1.659 

-2.558 

;  -1.643 

10.0 

-9.886 

|   -1.214 

-9.887 

;  -1.123 

-9.887 

-1.114 

50.0 

-49.98 

|   -1.120 

-49.98 

;  -1.030 

-49.98 

!  -1.021 

100. 

-99.99 

;   -1.101 

-99.99 

;  -1.020 

-99.99 

,  -1.011 

Table    6.        Non-Ideal  Pi-Section  Poles    (Figure    14). 
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Figure   V->.      ^on-Ideta   Fi^Section   Poles 
(See  Tnble   6) 


45 


synthesis.     But  solid  state  devices  are  extremely  vulnerable  to  the 
effects  of  radiation,    since  the  operation  of  most  of  these  devices  is 
determined  by  the  density  and  gradients  of  the  minority  carriers. 
This  minority  carrier  density  can  be  drastically  upset  by  radiation 
exciting  electrons  from  the  valence  band  to  the  conduction  band  of  the 
material.     When  a  transistor  in  a  circuit  such  as  figure  16  is  exposed 
to  a  high  energy  radiation  pulse  the  resultant  collector  current  wave- 
form is  similar  to  that  given  in  figure  17. 

The  integrated  circuit  transistor  operating  in  the  active  region 
can  be  considered  basically  as  one  forward-biased  and  one  reversed- 
biased  diode  connected  back-to-back.     When  irradiated  electron-hole 
pairs  are  generated  throughout  the  semiconductor  body  including  the 
substrate  (e.g.  ,    silicon).     In  the  transition  region  near  a  reverse- 
biased  p-n  junction  the  generated  minority  carriers  are  swept  across 
the  junction  by  the  electric  field  within  a  few  nanoseconds.     Carriers 
generated  outside  of  the  transition  region  will  raise  the  minority 
carrier  concentration  gradient  near  the  junction  generating  a  diffusion 
current.     These  two  components  of  current  constitute  the  primary 
photocurrent  in  a  p-n  junction  and;  thereby,   a  transistor.     The  general 
shape  of  the  primary  photocurrent  is  very  much  dependent  on  the 
minority  carrier  lifetime  as  indicated  in  figure  18.     In  an  npn  transistor 
the  generated  collector  region  minority  holes  flow  into  the  base  region. 
This  hole  current  is  amplified  by  the  current  gain  of  the  transistor 
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Figure  16.   A  dimple  Transistor  Circuit 
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Figure  17.   Transistor  Photocurrents  (induced) 
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Figure  18.   Relative  Shanes  of  the  Diffusion 

Component  of  Primary  Photocurrent. 
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circuit  generating  a  secondary  photocurrent  such  as  that  shown  in 
figure   17  . 

Radiation-induced  photocurrents  can  cause  errors  and  failures 
in  integrated  circuits.     If  the  dosage  is  large  enough  the  semiconductor 
material  can  suffer  permanent  displacements  effects  in  addition  to 
the  ionization  effects  which  are  usually  temporary    ["  27J      .     These 
effects  as  well  as  the  detailed  explanation  of  the  theory  and  relations 
between  the  radiation  dosage  (electron,    gamma,    etc.)  and  the  resultant 
photocurrent  magnitude  and  shape  we  leave  to  those  specializing  in 
such  matters.      Of  major  concern  here  is  the  .'fact  that  we  have 
"parasitic"  photocurrents  whose  effects  must  be  at  least  diminished 
if  not  eliminated.     One  method  is  circumvention,   which  in  its  most 
common  form  amounts  to  removing  power  from  critical  circuits 
during  the  radiation  pulse.     However,   in  many  applications  (e.g.  , 
counters,    guidance  mechanisms)  continuous  operation  is  mandatory. 
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Figure  19.       A  Simple  Logic  Gate 
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Dielectric  (e.g.  ,    silicon  dioxide)  isolation  is  a  very  effective 
means  of  removing  many  of  the  "evils"  of  radiation.     The  simple 
logic  gate  shown  in  figure  12  is  used  to  explain  the  effects  of  dielectric 
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Figure  20.        Logic  Gate  with  Parasitic  Elements 

isolation.     If  this  circuit  is  built  on  a  silicon  chip  by  a  diffused  process 
with  p-n  junction  isolation  (i.e.  ,   isolation  between  the  components 
and  the  substrate  is  obtained  by  reverse-biasing  the  p-n  junction 
boundaries),   the  parasitic  diodes  and  transistor  shown  in  figure  20 
will  generate  substantial  photocurrents  making  analysis  extremely 
difficult.     These  photocurrents  will  vary  widely  and  some  will  be 
quite  large.     Also  the  parasitic  transistor  forms  an  SCR  and  problems 
of  latch-up  are  encountered.     By  use  of  thin  film  resistors  and  di- 
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electric  isolation  that  is  relative  insensitive  to  radiation  (wide  for- 
bidden region)  these  parasitic  elements  and  their  leakage  currents 
can  be  eliminated  and  the  integrated  circuit  can  be  analysed  very- 
much  the  same  as  discrete  circuits  with,   of  course,   the  appropriate 
parameters.     The  resultant  equivalent  circuit  is  again  figure  19. 
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Figure  21.       Photocurrent  Generator  and  Compensation 

Diode  (dotted) . 

The  major  problem  is  now  the  transistor  photocurrent,   i       , 

PP 

as  shown  in  figure  21.     This  photocurrent  can  be  compensated  for  by- 
adding  a  diode  as  shown  (dotted).     If  thediode  junction  is  made  to  match 

the  collector-base  junction,   the  photocurrent,   i       . ,   of  the  compensating 

ppd 

diode  will  be  the  same  as  i        so  that  no  secondary  photocurrent  will 

PP 

result.     The  disturbance  due  to  the  transient  radiation  will  then  be 
kept  to  a  minimum. 


51 


No    C  ornpensa.t iop 


Compensation 

Photo  Current 


Figure  22.         Compensation  Definitions. 

By  defining  the  component  compensation  as 

I 


%  compensation    == 


N 


I 


(I       -     I      ) 

-£ EE_    (   100%) 


N 


I 


PP 


where  the  terms  are  defined  in  figure  22  it  is  possible  to  have  some 
measure  of  the  effectiveness  of  a  compensation  scheme  at  eliminating 
the  effects  of  the  secondary  photocurrents .     In  the  circuit  of  figure  21 
it  is  easily  seen  that  at  least  some  compensation  will  result  if  any 
compensating  diode  is  used  regardless  of  whether  or  not  the  collector- 
base  junction  is  identical  to  the  diode  junction.     Also  since  the  photo- 
current  of  a  p-n  junction  is  directly  proportional  to  the  junction  area 

[  25 J    the  ratios  of  i        and  i       ,  will  be  the  same  as  the  ratios  of  their 

PP  Ppd 

areas  with  the  amount  of  compensation  being  proportional  to  this 
ratio. 
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In  the  circuit  of  figure  21  the  amount  of  current  compensation 

will  be  determined  by  how  closely  the  junctions  can  be  made  to  match. 

However,    even  with  100%  compensation  this  configuration  suffers 

from  the  voltage  drop  of  i        through  the  collector  load  which,   in  itself, 

PP 

may  very  well  be  intolerable.     Also  as  pointed  out  by  Grannemann, 
et.   al.  ,    (J26J     the  amount  of  compensation  provided  is  very  sensitive 
to  a  close  match  of  the  compensating  diode  and  the  transistor  collector- 
base  junction.     This  sensitivity  is  due  primarily  to  the  amplification 
of  any  mismatch  by  the  gain  of  the  stage.     Utilizing  the  CIRCUS 
computer  program  [  12  &   1  3 J    an  analysis  of  figure  21  was  made 
using  the  parameters  of  the  1N903  and  2N404  as  "typical"  while  the 
compensation  diodes  were  special  formulations  to  match  (or  mismatch) 
the  collector  junction  of  the  2N404.      This  analysis  showed  that  even 
a  primary  photocurrent  of  0 .  3  milliamperes  peak  was  intolerable  in 
that  the  output  which  was  4.23  volts  with  zero  input  dropped  to  about 
zero  volts  and  stayed  there  for  about  four  microseconds  when  subjected 
to  a  10        rad  per  second  200  nanosecond  radiation  pulse  such  as  shown 
in  figure  17 . 

By  attaching  an  additional  identical  compensation  diode  across 
the  load  resistor  as  shown  in  figure  23  we  were  able  to  achieve  sub- 
stantially better  compensation  in  that  V  remained  more  nearly 

out 

constant.     The  sensitivity  of  the  output  to  a  mismatch  of  this  load 

diode  is  much  less  than  the  base  as  there  is  no  gain  following  it. 

But  as  noted  above  imperfect  matching  of  the  base  diode  to  the  transistor 
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Figure  23.       Base  and  Load  Diode  Compensation 

is  magnified.     Next  over-compensating  with  the  load  diode  was  tried. 

This  led  to  an  "increasing"  output  voltage  (with  zero  input)  during  a 

pulse.     This  factor  was  immediately  seized  on  as  a  means  to  decrease 

the  sensitivity  of  the  circuit  to  the  amplified  mismatch  of  the  base 

diode.     The  amplified  mismatch  lowers  V         while  the  over  compensated 

out 

load  diode  raised  V        .     This  over- compensation  can  be  effectively 

out 

accomplished  by  increasing  junction  area  of  the  load  diode.  By 
utilizing  a  factor  of  2  and  our  values  a  compensation  of  80%  was 
attainable  for  a  base  diode  tolerance  of  20%. 

Next  the  circuit  of  figure  24  was  attacked.     Initially  only  base 
diodes  on  Tl  and  T2  were  utilized  and  theoretical  100%  compensation 
was  attainable  for  each  transistor  if  the  diodes  and  transistors  were 
exact  matches.     However,    sensitivity  to  mismatch  worsened  as  Tl 
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Figure  2i|.   The  Gate  Circuit. 
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Figure   25.      Condensations    for   the   Gate   Circuit. 
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mismatch  is  amplified  by  T2  when  T2  is  turned  on  by  its  own  mis- 
match.    As  a  tolerance  of  less  than  1-5%  is  not  feasible  in  integrated 
circuits  with  the  present  technology,   this  approach  was  abandoned. 

As  T2  forms  part  of  the  emitter  circuit  of  Tl  attempts  to  com- 
pensate solely  T2  met  with  very  little  success.     By  attaching  base 
diodes  to  both  Tl  and  T2  along  with  load  diodes  across  the  1  .75  K  and 
the  6.0  Kohm  resistors  as  shown  in  figure  25  adequate  compensation 
was  effected.     As  mismatch  at  Tl  would  cause  a  secondary  photo- 
current  into  the  base  of  T2  the  base  diode  of  T2  should  be  slightly 
greater  in  junction  area  than  Tl  base  diode  to  provide  maximum  area 
matching  tolerance.     With  the  factor  of  two  between  base  and  load 
diode  junction  areas  a  compensation  of  80%  was  attainable  with  areas 
tolerance  of  20%. 

A  theoretical  development  of  an  effective  and  economical 
means  of  radiation  compensation  in  integrated  circuits  (also  discrete 
but  not  so  economical)  has  been  demonstrated.     An  iterative  transient 
analysis  is  required,   thereby,   virtually  demanding  the  use  of  a 
computerized  analysis.     The  costs  associated  with  constructing, 
testing,   and  evaluation  actual  devices  without  the  use  of  such  an 
analysis  would  be  prohibitive  while  with  the  analysis  at  least  a  reason- 
able first  cut  can  be  made. 
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APPENDIX  A 
CLASS  "C"  AMPLIFIER 

The  following  program  selects  appropriate  parameter  values 
to  optimize  a  class  "C"  circuit.     As  rather  complete  comments  are 
contained  in  the  listing  a  detailed  explanation  is  not  included  at  this 
point,   but  two  points  should  be  brought  out.     The  optimizing  function 
is  calculated  in  subroutine  OPT  and,   therefore,   it  is  easily  modified 
to  select  or  weigh  some  other  parameters.     Similiarly  KCURR 
describes  a  tube  model  that  can  be  quite  easily  changed  to  some 
other  model  (e.g.  ,   a  transistor  model). 
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MAP: 


APPENDIX  B 


FILEMAP  AND  FORMAT 


CIRCUS  DISTRIBUTION  TAPE 


File  1  CIRCUS  OBJECT  DECK:    From  FORTRAN  H  OPT 

0  COMPILER  and  release  9  assembler.     4000- 
5000  cards.     May  be  loaded  directly  from  tape. 

File  2  CIRCUS  MODEL  LIBRARY  UPDATE  PROGRAM 

OBJECT  DECK  (complete)  approximately  1000 
cards. 

File  3  CIRCUS  FORTRAN  SOURCE  DECK.     Approximately 

6000  cards. 

File  4-12  CIRCUS  ASSEMBLER  LANGUAGE  ROUTINES: 

Approximately  500  cards  total. 

File  13  CIRCUS  MODEL  LIBRARY  UPDATE  PROGRAM 

FORTRAN  SOURCE  DECK,    exclusive  of  those 
routines  common  to  CIRCUS:     See  CIRCUS  PRO- 
GRAM MANUAL  (No.    346-2)  less  then  500  cards. 

File  14  CIRCUS  SAMPLE  DATA.     Approximately  500  cards. 


FORMAT: 


The  14  files  fill  approximately  1000  feet  of  tape. 

All  files  on  all  tapes  are  80  character  card  images:     80  byte 

records  unblocked. 
All  tapes  are  odd  parity,   density  800  bits/inch. 
Note  that  the  first  two  files  only  are  in  object  form,   all  others 

are  in  EBCDIC. 
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NOTES  ON  CIRCUS/ 360: 

In  order  to  CIRCUS/ 360  as  supplied,   the  user  must  provide 
timing  routines  which  are  called  by  the  program.     These  routines  are 
a  part  of  Lincoln  Laboratory's  system  library  and  therefore  are  not 
a  part  of  CIRCUS.     Dummy  routines  may  be  written,   but  if  different 
names  are  used,    or  if  similar  routines  exist  in  the  user's  system  but 
with  different  names,   the  calling  statements  must  be  changed  in 
CIRCUS.     Calls  appear  in  KLOCK  and  OUTPUT  in  several  places  in 
these  routines.     The  routines  are: 

TIMDAT    (IYR,   IDAY) 

TIMHR    (HOURS) 

TIMHMS    (IHR,   IMIN,  ISEC) 

where  I  implies  fixed  point  and 

H  implies  floating  point  numbers 

These  routines  use  the  system  clock  and  times  as  implied  in  their 
arguments  . 

If  the  system  on  which  CIRCUS/  360  is  to  be  used  has  a  core 
capacity  greater  than  256  K  bytes,   the  capacity  of  CIRCUS  may  be 
increased  to  handle  circuits  of  greater  than  70  equations  by  adding 

C0MM0N  Al(50000)  CIR25715 

and  changing  CIR25750  to  read 

LC0RE    =     50000  CIR25750 

in  SUBROUTINE  MAPPER. 
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In  order  that  the  supplied  version  fit  into  a  256  K  byte  core, 
the  system's  load  point  must  be  B000  or  less. 

CIRCUS/ 360  uses  the  system  clock  via  the  assembler  language 
routine  SETLMT  to  compare  execution  time  against  a  specified  limit. 
If  no  limit  is   specified,   the  program  assumes  a  ten  minute  limit, 
after  which  CHARET  is  called  and  execution  terminated.     To  adjust 
this  limit,   the  first  card  in  the  data  deck  should  read,   beginning  in 
column  one , 

LIMITS,   X 
where  X  is  the  desired  limit  in  minutes,   floating  point  (as  in  LIMITS, 
8.5  for  8  min.    30  sec.   maximum  execution). 

CIRCUS/ 360  now  includes  a  printer  plotter  routine  which  is 
called  by  the  PLOT  statement  in  the  input  list  (see  user's  manual). 
Plots  are  51  by  120  characters;  time  interval  is  the  smallest  in  the 
INTERVALS  statement;  if  greater  than  120  time  points  at  this  interval 
exist,   the  plot  will  be  continued  on  succeeding  pages.     One  plot  vari- 
able appears  per  plot. 

KNOWN  BUGS  IN  CIRCUS/  360 

1.  Occasionally  a  circuit  time  and  associated  outputs  may  be 

printed  twice.     This  occurs  at  times  when  something  is 
specified  in  the  input  list  as  changing.     The  first  of  these 
printed  times  represents  a  value  slightly  less  than  that 
printed  by  an  amount  determined  by  machine  accuracy,   and 
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rounded  off  when  printed.     The  second  value  is  the  exact  time 
as  specified  in  the  input  list. 
2.  When  specifying  a  logical  unit  for  a  model  library,    save  tape, 

etc.  ,   a  user  is  supposed  to  be  able  to  specify  the  logical 
unit  number  (as  in  SAVE,    15)  or  use  a  letter/  number 
equivalent  (as  in  SAVE,   E)  where  the  letter/number  cor- 
respondence is  specified  in  LINK2  (DATA  MNSAVE/1HA, 
1HB,1HC,  JHD,  1HE,  1HF,  1  ,2,3,4,  15,  16/  CIR20600). 
As  it  now  stands,   the  number  designation  will  not  be 
recognized,   while  the  letter  designation  works.     If  a  number 
other  than  those  in  the  standard  letter/ number  correspondence 
is  desired,   the  aforementioned  DATA  statement  may  be 
changed  to  provide  a  more  useful  correspondence. 
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APPENDIX  C 

THE  CODING  OF  POTTLE 
LINEAR  NETWORK  ANALYSIS  PROGRAM 

Given  an  arbitrary ,    consistent,   linear,   fixed  electrical 
network  (with  zero  initial  conditions),    the  POTTLE  program  cal- 
culates: 

(a).     A  prescribed  but  arbitrary  transfer  function  t(s)  or, 

more  generally,   an  arbitrary  n-port  matrix  t..(s). 
(b).     The  poles  and  zeros  of  t(s)  or  t..(s). 
(c).     The  frequency  responses  associated  with  t..(s)  - 

prints  and  printer-plots  magnitude,   phase,   and  delay. 
(d).     A  set  of  state  equations  for  the  network. 

In  addition,   any  set  of  element  values  may  be  stepped  up  to 
six  times  and  the  above  analysis  repeated  at  each  step. 
The  network  size  limitations  are: 
(a).     Maximum  number  of  elements  -   60 
(b).     Maximum  number  of  LC  elements  -   32 
(c).     Maximum  number  of  nodes  -   27 
(d).     Number  of  controlled  sources  plus  number  of  inputs 

must  be  less  than  or  equal  to  32. 
(e).     Number  of  controlled  sources  plus  number  of  outputs 

plus  number  of  resistors  must  be  less  than  or  equal  to  64, 
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(f).     Three  time  (number  of  LC)  plus  two  times  (number  of 
inputs)  plus  number  of  branches  minus  number  of 
nodes  must  be  less  than  or  equal  to  128. 

(g).     Maximum  number  of  stepped  element  values  -  88 


The  input  data  format  is: 

(a).     A  title  card  -   columns   1-80  will  be  printed  on  the 

output  sheet. 
(b) .     Cards  describing  the  branches  of  the  network: 
Columns  Des  cription 


1   -  4 


6-10 


11 
14 


12 
15 


Branch  name  (left  justified)  beginning 
with: 

V  -  voltage  source 
I    -   current  source 
R,    L,    or  C  -   resistor,   inductor, 
or  capacitor,    respectfully 

Name  of  controlling  branch  (left  justified) 
if  col  1-4  contains  a  controlled 
source;  otherwise  leave  blank 

Positive  node  (right  justified) 

Negative  node  (right  justified) 


Column 
11-12       O- 


Column 
-O       14-15 


Number  nodes  consecutively  from  zero 
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Columns  Description 

18  If  branch  variable  is  to  be  considered 

an  output  quantity,    use  symbols  V 
or  I  as  appropriate. 

20  Type  of  control  (V  or  I)  if  controlled 

source;  otherwise  blank 

21     -     30  Value  of  element  or  strength  of  con- 

trolled source;  independent  source 
leave  blank.  Must  have  a  decimal 
point  (F10.0  format) 

31     -     40  Values  of  stepped  elements  or 

strengths  of  stepped  controlled 

41     -     50  sources;  all  stepped  elements 

must  be  stepped  the  same  number 
of  times;   stepped  element  val- 
ues may  not  be  continued  on 

71     -     80  succeeding  cards.     Must  have 

decimal  point  (F10.0  format) 


(c).     A  card  denoting  the  types  of  responses  desired.     This 

card  must  be  present. 

columns   11     -     12:      01  if  frequency  response  is 

desired 

00  or  blank  otherwise 

(d).     Two  cards  to  describe  frequency  response  desired;  these 

are  omitted  if  no  frequency  response  is  desired. 

Card  1: 

column  1:         1  if  linear  frequency  scale 
2  if  log  frequency  scale 

column  3:         1  if  magnitude  calculated  in  db 

2  if  magnitude  in  ohms 

3  if  magnitude  in  mhos 
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column  linear  frequency  scale 

total  number  of  frequencies 

5     -     7: 

logrithmic  frequency  scale 

number  of  frequencies  per 
decade 

number  is  right  justified 
(13  format) 

Card  2 
columns      1   -    10:     Lowest  frequency  (Hz) 
columns   11   -   20:     Highest  frequency  (Hz) 

(e).     If  more  than  one  network  is  to  be  analyzed,   a  set  of  these 
cards  is  required  for  each  additional  network. 

Additional  comments: 

(a).     The  method  of  analysis  used  is  sensitive  to  scaling, 

producing  overflows  if  improperly  scaled.     Values  for 
L's  and  C's  should  be  within  an  order  of  magnitude  of 
unity  for  large  networks. 

(b) .      Zero  values  for  inductors  and  capacitors  may  be  used 
to  establish  a  dummy  branch  current  or  node-pair 
voltage;  this  is  useful  in  controlling  a  controlled  source 
or  defining  an  output  variable. 
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APPENDIX  D 
THE  POTTOL  TRANSLATOR 

PURPOSE 

The  primary  purpose  of  the  Translator  is  to  convert  an 
engineer-oriented  free-format  circuit  description  into  the  standard 
input  format  required  by  C.   Pottle's  Linear  Network  Analysis 
Program.     Prior  exposure  to  the  POTTLE  program  (See  Appendix 
C)  is  assumed. 

FEATURES 

The  major  features  of  the  Translator  are: 

(1).  Free-format  of  input  data 

(2).  Minimal  changes  to  the  original  program 

(3).  Simple  and  effective  overlay  (if  needed) 

(4).  Additional  diagnostics 

(5).  Exact  reprint  of  input  data 

(6).  Continuation  cards 

(7).  Comment  cards 

(8).  Batch  processing 

OVERLAY 

The  particular  form  in  which  the  Translator  has  been  con- 
structed lends  itself  to  a  very  natural  and  fast  overlay  structure. 
Figure  D- 1    shows  the  recommended  overlay  structure.     This  particular 
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MAIN  PROGRAM 
SUBROUTINES 
INPUT 
DECODE 


TRANSLATOR 
LANPRO   FREPRO 
TNUM     OUTPUT 
PACK     TAPE 
EDIT      SHIFT 


Place  all  other  routines 
in  this  branch 


SHIFT  is  in  Assembly  language 


Figure  D- 1  .        Overlay  Tree  Structure. 

structure  loads  the  Translator  along  with  the  root  segment,   processes 
the  input  data,   loads  the  rest  of  the  program,    executes  the  program 
and  stops  without  having  to  repetitiously  load  and  reload  the  same 
overlay  branch. 

POTTLE  PROGRAM  CHANGES 


Since  the  logical  unit  assignments  used  in  the  original  program 
were  not  the  same  as  used  locally,   it  was  judged  worthwhile  to  change 
all  logical  unit  references  to  variables  and  assign  them  at  one  central 
location.     This  is  presently  accomplished  in  the  main  program  of  the 
program.     These  assignments  are  transferred  via  labeled  COMMON 
to  the  various  subroutines.     The  present  assignments  and  their 
correlation  with  the  original  assignments  are  shown  in  figure  D-2. 
These  assignments  may  be  changed  permanently  by  changing  the 
appropriate  cards  in  the  main  program.     Additionally  a  temporary 
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Variable 

Unit 

Use 

Original 

INTAPE 

5 

User's  input  data 

1 

IOUTP 

6 

Copy  of  input  and  errors 

2 

STATP 

6 

State  equations 

3 

TRANTP 

6 

Transfer  functions,   poles/zeros 

4 

FREQTP 

6 

Frequency  response 

5 

NOT 

E:      Um 

it  6  is  the  printer  locally 

Figure  D-2.        I/O  Unit  Assignments. 

change  of  any  or  all  of  units  may  be  invoked  for  a  single  run  through 
the  use  of  a  TAPES  card  as  shown  in  figure  D-3. 


column  7  or  later 

TAPES,  STATP    =     2,  TRANTP    =  3 


Figure  D-3.       Temporary  Unit  Assignments 


ADDITIONS 


The  additional  features  added  to  the  allowable  input  include: 
(1).     Continuation  Cards:     Any  number  of  continuations  are 
allowed  as  long  as  the  previous  card's  last  character 
is  a  comma.     There  is  a  limit  of  154  non-blank  charactors 
after  column  6  of  the  first  card. 
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(2).     Debugging:     A  limited  debugging  facility  is  provided  in 
that  a  card  containing  the  word  TRACE  (starting  in 
column  7  or  later)  will  cause  a  diagnostic  print  as  control 
enters  and  departs  each  subroutine.     It  also  provides  a 
print  out  of  the  NHEX  variables  in  Al  format. 

(3).   Batch  Processing:      Batch  processing  is  implemented 
by  requiring  that  each  set  of  circuit  description  data 
end  with  an  EXECUTE  card  and  the  batch  run  or  job 
end  with  an  END  card. 

(4).  Comment  Cards:  Except  for  reproduction  in  the  print 
out,   any  card  with  an  asterick  in  column  one  is  ignored. 

(5).  Number  Format:  All  numbers  on  the  input  data  cards 
may  be  either  F,  E,  or  integer  format  with  or  without 
decimals.     Recognizable  number  forms  are: 

-XX.XE-YY  XX  .XXEYY 

XXXE+YY  X.X  +  .XX 

If  the  decimal  is  missing  it  is  assumed  to  follow  the 
last  X. 

(6).   A  print  out  of  the  translated  program  is  provided. 
INPUT  DESCRIPTION 


The  standard  R,  L,  C,V,   or  I  branch  is  coded  by  placing  a 
unique  name  for  the  branch  in  columns  1-5  (note:     This  is  the  only 
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type  of  card  allowed  to  have  anything  in  the  first  five  columns).     If 
five  characters  are  used  in  the  name  the  fifth  character  will  be 
ignored;  however,   the  name  need  not  start  in  column  one.     All  other 
information  is  placed  after  column  six  (including  continued  cards). 
The  node  numbers  are  delimited  by  parens  and  separated  by  a  comma, 
The  order  of  the  node  numbers  implies  the  positive  direction  of 
current  flow.     Spaces  may  be  used  freely  to  enhance  readability. 
Typical  examples  are  shown  in  figure  D-4.     In  the  case  of  stepped 
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Figure  D-4.       Input  Data  Examples. 

elements  the  initial  value  (or  only  if  not  stepped)  is  placed  first  with 
sequential  stepped  parameter  values  separated  by  commas  following. 
All  stepped  parameters  must  be  stepped  the  same  number  of  times. 
An  output  of  the  branch  voltage  or  current  is  indicated  by  adding 
the  expression  OUTPUT  =  V  or  OUTPUT  =  I,    respectfully,   to  the 
branch  description. 

Descriptions  of  controlled  sources  have  one  additional  restric- 
tion.    The  controlling  branch  name  must  appear  BEFORE  the  type 
of  control.     This  is  illustrated  in  figure  D-5. 
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Figure  D-5.       Description  of  Controlled  Sources 

FREQUENCY  RESPONSE 

The  inclusion  of  a  card  of  the  form  shown  in  figure  D-6  will 
generate  the  proper  commands  to  obtain  the  indicated  frequency- 
calculations  .     This  card  will  set  the  variables  via  labelled  common 
required  in  subroutine  TPRINT.     It  is  not  necessary  to  include  the 
words  LOG  and/or  DB  if  these  options  are  desired. 

COMMANDS 


The  EXECUTE  card  causes  an  analysis  to  begin.  It  also 
generates  the  blank  card  required  by  the  POTTLE  program  if  no 
frequency  card  is  included. 


FREQ  =  nl     (n2)     n3,       LINEAR  or  LOG,   0HMS  or  MH0S  or  DB 
nl-min  freq  (Hz  Defaults 

n2-no.   freq/ decade  (log)  FREQ 

no.   freq  (linear  freq  scale 

n3-max  freq  (Hz)  units 


none 
LOG 
DB 


Figure  D-6.         The  Frequency  Card, 
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An  END  card  terminates  the  particular  run.     Lack  of  an  END 
card  does  not  adversely  affect  the  execution  of  the  program  (if  all 
other  data  is  correct). 

COMMENTS 

The  Translator  is  intimately  tied  to  the  POTTLE  program. 
Therefore  a  proper  appreciation  of  the  Translator  requires  at  least 
a  passing  acquaintance  with  the  use  of  the  original  POTTLE  program. 
To  this  end  Appendix  C  contains  a  description  of  the  original  program 
and  the  CAD  LIB  contains  the  original  program. 

In  order  that  POTTLE  changes  could  be  kept  to  an  absolute 
minimum  the  Translator  operates  on  only  one  card  (and  its  con- 
tinuations) at  a  time,   i.e.  ,   the  order  of  input  is  the  same  as  the 
original  POTTLE  program  except  that  the  label  card  is  no  longer  used 
as  unlimited  comment  cards  are  available  for  expositor  comments. 


APPENDIX  E 
CALAHAN    -     A  USER'S  GUIDE 

CIRCUIT  LIMITATIONS 

(1)  Complexity 
(a) .     30  nodes 

(b).      100  passive  elements 
(c).     20  active  elements 

(2)  Voltage  controlled  current  source  is  the  ONLY 
allowable  active  element. 

(3)  Storage  -  approximately  90  K  bytes 

DEFINITIONS 

(1)    Network  -  the  circuit  to  be  analyzed  when  considered 
as  a  two-port  black  box  (See  figure  E- 1  below). 
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Figure  E-l.       The  Calahan  Network. 


(2)    Network  function  -  A  functional  relationship  between  the 
input  voltage  (V   ),   the  output  voltage  (V   ),   the  input 
current  (I   ) ,   and/ or  the  output  current  (I   )  where  these 
terms  are  positive  as  indicated  in  figure  E-l. 
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(3)    KEY  1   -  A  computer  flag  whose  value  indicates  a  particular 
network  function  (See  Table  E-l  below). 


KEY  1 

Network  Function 

Symbo 

lie 

1 

Voltage  Transfer  Function 

V, 

I2  =  0 

2 

Open  Circuit  Driving 
Point  Impedance 

v, 

J2  =0 

3 

Open  Circuit  Transfer 
Impedance 

I, 

!*=* 

4 

Short  Circuit  Driving 
Point  Impedance 

Vz=* 

5 

Short  Circuit  Transfer 
Admittance 

v, 

V2  -0 

6 

Current  Transfer  Function 

I* 

V2  =0 

Table  E-  1 .       Definitions  of  KEY  1 . 

(4)    KEY  2  -  A  computer  flag  whose  value  indicates  a  particular 
set  of  desired  outputs  (See  table  E-2). 


INPUTS 


(1)  Network  elements  or  network  function  (Voltage  or  current 
transfer  function  only. 

(2)  Input  and  output  node  numbers 

(3)  Type  of  network  function  desired  as  output 
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KEY  2 

Output(s)  is  (are) 

of  the  Network  Function 

1 

Coefficients,   poles, 

and  zeros 

2 

Coefficients,   poles, 

2eros,   and  freq  response 

3 

Coefficients,   poles, 

2eros,   and  time  response 

4 

Coefficients,   poles, 

zeros,   freq  response,   and 

time  response 

5 

Symbolic  form  of  the  network  function 

Table  E-2.       Definitions  of  KEY  2. 

OUTPUTS  (See  Table  E-2) 

(1)  Coefficients  of  requested  network  function 

(2)  Poles  and  zeros  of  the  requested  network  function 

(3)  Frequency  response  (both  tables  and  printer  plots  of 
magnitude,   phase  and  delay) 

(4)  Time  response  to  an  arbitrary  input  (combination  step 
in  CADLIB  version).     Tables  and  printer  plots  are 
included. 

(5)  Repeated  outputs  of  1-4  with  some  element  value  as  a 
variable  parameter 

(6)  Designated  network  function  in  symbolic  form 

INPUT  CODING 

The  initial  step  in  coding  for  CALAHAN  is  to  number  all  nodes 
sequentially  from  one.     A  deck  of  cards  is  then  coded  as  indicated  on 
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the  flowchart  of  figure  E-2.     Two  points  need  clarification.     The  first 
card  is  a  title  card  and  is  ignored  except  for  titling  the  output.     The 
second  card  is  a  dual  purpose  card.     If  the  input  is  in  network 
element  form,    card  #2  contains  a  positive  number  equal  to  number 
of  circuit  descriptions  between  title  cards.     As  each  circuit  or  circuit 
variation  will  normally  be  different  the  usual  value  will  be  one.     If 
the  input  is  to  be  a  transfer  function  the  value  used  will  be  a  negative 
KEY  2. 

DISCUSSION 


General  Precautions 

(1).     Poles,   zeros,   and  time  response  may  be  impossible  to 

find  for  no  predictable  reason. 
(2).     Some  bandpass  networks  may  yield  incorrect  responses 

due  to  single  precision  limitations. 
(3).     The  nodes  must  be  numbered  1  through  N,   the  number 

of  nodes;  otherwise,   there  are  no  constraints  on  the 

manner  of  numbering. 
(4).     The  network  must  be  connected  (i.e.  ,   a  path  must  exist 

between  any  two  nodes).     Thus,   for  a  transformer,   the 

common  connection  between  windings  must  be  made, 

unless  a  path  connecting  the  windings  already  exists 

through  external  circuitry. 
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(5).     The  RLC  elements  must  be  listed  before  the  controlled 

sources,   although  in  any  order. 
(6).     Two  or  more  RLC  elements  across  the  same  nodes  must 

have  nodes  listed  in  the  same  order, 
(7).     Negative  RLC  elements  are  permitted. 
(8).     A  zero  value  for  a  resistance  or  inductance  is  never 

permitted. 
(9).     If  repeated  analyses  are  made  with  some  element  value 

as  parameter,   a    zero  value  is  never  permitted  for  a 

resistance  or  inductance,   and  is  permitted  only  as  a 

final  value  for  a  capacitance  or  controlled  source. 
(10).    When  a  logrithmic  frequency  scale  is  used,    recall  that 

zero  frequency  is  not  a  valid  lowest  frequency. 

(11).    When  choosing  a  highest  frequency,  beware  of  possible 

overflow.     If  a     is  the  coefficient  of  the  highest  power 
n 

(n)  of  s  in  the  denominator,   then  overflow  may  occur  if 

O  no 

(2  31  f    a   )         is  greater  than  (10) 

(1Z).     The  time  response  cannot  be  calculated  for  a  network 
function  with  a  numerator  of  higher  degree  than  the 
denominator . 

(13).    The  topological  program  can  be  fooled  by  certain  de- 
generate cases  which  do  not  warrant  safeguarding. 
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